clear;
[p,s,d] = Msi3dReadData('27010000_0820_X_P45_10T_P2 _3D.txt');
[phi2,theta2]=meshgrid(p*pi/180,s*pi/180);

[x,y,z]=sph2cart(phi2,theta2,d);
clear p d s;
%subplot(2,1,1);
mesh(x,y,z);
%subplot(2,1,2);
%surf(x,y,z);